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Abstract. We present near-infrared J band photometric observations of the intermedi- 
ate polar WX Pyx. The frequency analysis indicates the presence of a period at 1559.2 
+ 0.2 seconds which is attributed to the spin period of the white dwarf. The spin period 
inferred from the infrared data closely matches that determined from X-ray and optical 
observations. WX Pyx is a system whose orbital period has not been measured directly 
and which is not too well constrained. From the IR observations, a likely peak at 5.30 
+ 0.02 hour is seen in the power spectrum of the object. It is suggested that this corres- 
ponds to the orbital period of the system. In case this is indeed the true orbital period, 
some of the system physical parameters may be estimated. Our analysis indicates that 
the secondary star is of M2 spectral type and the distance to the object is 1.53 kpc. An 
upper limit of 30° for the angle of inclination of the system is suggested. The mass 
transfer rate and the magnetic moment of the white dwarf are estimated to be (0.95 - 
1.6)xlO" 9 M yr 1 and (1.9 - 2.4)xl0 33 G cm 3 respectively. 
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1. Introduction 

Intermediate polars (IPs) are a sub-class of magnetic cataclysmic variable stars (mCVs). An IP 
comprises of a semi-detached binary system consisting of a white dwarf primary star and a late- 
type Roche-lobe filled secondary s tar; the white dw arf is believed to have a moderate magnetic 



field strength of B ~ 1 - 10 MG dPatterson 19941) . The rotation period of the white dwarf is 



not synchronized with the binary orbital period. The mass transferred from the secondary to the 
primary, by the process of Roche lobe overflow, eventually accretes onto or near the magnetic 
poles of the white dwarf through the magnetic field lines. The shock heated plasma produces 
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X-rays by the cooling o f electrons near the white dwarf's surf ace through free-free interactions 
( Evans & Hellierll200 7i). Recently. iBednarek & Pabich ( 2010b have proposed that y rays could 
also be produced when the accelerated hadrons are convected to the white d warf surface a nd 
interact with the dense matter. A comprehensive review of IPs may be found in lWarnerl (11995b . 



WX Pyx was identified as IE 083 0.9-2238 in the Galactic plane survey of Einstein X-ray 
observatory by iHertz & Grindlav ( 1984 1. The blue excess seen in the object and also the nature 
of the optical spectra, which displayed emission lines including Hydrogen Balmer lines, Hei 
and Hen - lines which are typically seen in the spectra of CVs - suggested the object as a CV 
candidate (Hertz et al. 1990). The strong presence of the high-excitation Hen 4686 A line in the 



optic al spectrum also suggested the possibility o f a magnetic nature for the system (Hert z et al. 



opti 

1990). From an analysis of the optical light curve lO'Donoghue ( 19961) found a stable period of 
26 minutes which was attributed to the spin period of the white dwarf and confirmed the IP nature 
of the WX Pyx system. They suggested a value of of ~ 6 to 9 hours for the orbital period, based 
on spectroscopic and photometric data, but stressed that this result was not robust and should be 
viewed with caution. Though there are no direct X-ray observations of WX Pyx available till 
date, the object was accidentally detected 10.5 arc -minutes off-axis during a pointed observation 
of NGC 2613 by XMM -Newton. ISchlegell d2005l) analyzed data related to this serendipitous X- 
ray detecti on and found a spin pe riod of 1557.3 s which matched well with the optical spin-period 
estimate of lO'Donoghud (119961) . An orbital period of ~5.54 hours was inferred indirectly from 
the separation of spin-orbital side band frequencies in the power spectrum. 



The present study is partially motivated by the fact that there are no infrared studies of WX 
Pyx. More importantly, we had hoped to determine more robustly the orbital period of the system 
for which no direct observational evidence is available. Such a determination would allow many 
of the system parameters and mode of accretion of WX Pyx to be better estimated. In this study 
we present photometric observations of near-IR J band and undertake a time series analysis to 
estimate the orbital and spin period of the object. 



2. Observations and data analysis 



Near infrared J band photometry of WX Pyx was done using the Mt.Abu 1.2-m telescope. The 
object was observed for a total duration of 22 hours spanning six nights from December 2007 
to February 2010. The observations were taken using the Near-Infrared Imager/Spectrometer 
which uses a 256x256 HgCdTe NICMOS3 array with a FOV of 2x2 arc-min 2 . The telescope 
was dithered a t 5 diff erent positions during observations, as is customary for IR observations (for 
e.g. Das et al. ( 2008b ). to produce median sky-cum-dark frames for individual images. Several 
selected field stars were always kept in the field of view while dithering to enable differential 
photometry to be done. Frames of smaller durations were co-added when single frames of lar- 
ger duration were not available. Flat fielding, using twilight flats, was done to ensure a proper 
pixel response of the detector. An appropriate median sky-cum-dark image was generated and 
subtracted from the object frames. Instrumental magnitudes, obtained from these sky-subtracted 
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Table 1. Log of the photometric observations. 



Date of 


Reduced 


Filter 


Exp. 


Time 


Spin 


observation 


JD 




time 


cover- 


cycles 


(dd/mm/yyyy) 






(sec) 


age (hr) 


covered 


18/12/2007 


54452 


J 


180 


4.5 


10 


19/12/2007 


54453 


J 


180 


3.9 


9 


18/01/2008 


54483 


J 


60 


3.3 


7.5 


19/03/2009 


54909 


J 


60 


4.1 


9.5 


19/03/2009 


54910 


J 


60 


4.0 


9.5 


22/02/2010 


55250 


J 


180 


2.2 


5 



Note: Reduced JD = JD - 2400000 



frames using aperture photometry, were finally flux calibrated by comparing with several field 
stars whose lightcurves had previously been checked for photometric stability. 

The routines used for data reduction are based on the Interactive Data Language (IDL) pack- 
age. Specifically, APER procedure from GSFC IDL astronomy user's library was used for syn- 
thetic aperture photometry. The log of the observations is given in Table 1, which gives the date 
of observation, filter information, exposure time, total observation coverage and number of spin 
cycles covered. 

3. Results and discussion 
3.1 Light curves and frequency analysis 

The light curves of individual nights are shown in Figure 1 along with their corresponding peri- 
odograms. A clear modulation at a period close to 26 minutes can be visually seen in all the 
light curves. Equivalently, analysis yields the presence of strong peak near 0.64 mHz in all the 
periodograms. The mean / magnitude of the program star, peak frequency and semi-amplitude of 
the light curve on individual nights are presented in Table 2. The peak frequency varies between 
0.6367 mHz to 0.6439 mHz in periodograms of individual runs. This variation of 0.007 mHz, 
however, is much below the error level estimated from half width at half maximum of the peak 
which is 0.02 mHz. 

Rapid but weak oscillations at frequencies other than that of the principal component at 0.64 
mHz were also observed in certain light curves. Notable among these is the one near 1.2 mHz 
which is the first harmonic of the spin period. In order to check the coherence of these weak 
oscillations we plotted the periodograms up to the Nyquist frequency for all the light curves of 
individual nights. We do not find any coherent frequency above 0.64 mHz. We thus conclude that 
these weak oscillations may be due to flickering in the light curve, a feature which is a common 
characteristic of CVs. 
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Table 2. Frequency analysis of daily lightcurves. 
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Figure 1. Light curves of all the nights listed in Table 1 are shown in the panels on the left. The abscissa 
is HJD - HJD0 where HJD0 is arbitrary and given in the lower left for each day. Solid line shows the best 
fit sinusoid plus a polynomial fit representing the long term variation. The Lomb-Scargle periodograms of 
individual runs are shown in the right hand panels. It may be noted that the periodograms extend up to the 
Nyquist frequency. 
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Figure 2. The Lomb-Scargle periodogram of all the data is shown. The dotted lines indicate orbital and 
spin frequencies. The dashed lines show the first harmonic of the orbital and spin periods whereas the dot- 
dashed lines represents the spin-orbital sidebands. The window function is plotted as an inset in the top of 
the figure. 

The periodogram of the combined data plotted till the Nyquist frequency is shown in Figure 
2. A strong peak can be seen near 0.64 mHz. Another noteworthy peak is at the low frequency end 
of the periodogram near 50 /iHz which is discussed subsequently. One can also see minor peaks 
corresponding to 0.74 mHz and 1.28 mHz. But we consider these latter peaks to be statistically 
insignificant because the associated power is at a level of only 2 sigma above the background. 
Temporal gaps of various lengths in the observed data set cause the various aliasing peaks to be 
present in the periodogram of Figure 2. One day aliases are dominant in particular. This can be 
clearly seen in a plot of the window function which is presented as an inset in Figure 2. The 
periodogram between frequencies 0.4 and 0.9 mHz is shown in Figure 3 which is magnified part 
of the region containing peak near 0.64 mHz frequency. The highest power is at a frequency of 
0.64134 ± 0.00008 mHz which corresponds to a period of 1559.2 + 0.2 seconds. It may be noted 
that the formal error of the frequency (or period) is estimated by half-width at half-maximum of 
the highest peak in the peridogram which is the conventional way to determine the error. 

It is possible to generate the pulse profile of the spin modulation from the entire data set. As a 
first step, a best fit polynomial of either degree one or two was subtracted from data of individual 
night's to remove the long term variation from the light curve. The data were then folded at a 
period 1559.2 second. The phase was divided into 25 bins to produce the spin pulse profile which 
is presented in Figure 4. The error in each bin is inferred from the standard deviation of the data 
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Figure 3. The LS periodogram of all the data near the spin frequency, which is marked with a vertical dash 
at 0.64134 mHz, is shown. 



points within the bin. From our analysis we find that the time of maximum of the oscillation is at 
RJD 54322.03408 + 0.00005 d. 

In Figure 2, another significant peak is seen at a much lower frequency near 50 //Hz. A Lomb- 
Scargle periodogram produced around this frequency, in the range 0.02 //Hz to 200 yuHz, is shown 
in Figure 5. In this periodogram the maximum power is present at a frequency of 52.43075 /Aiz 
corresponding to a period of ~5.3 hours. However, various alias peaks are strong in this spectrum 
and the 1-day alias in particular, at a frequency of 40.45529 //Hz or alternatively corresponding to 
a period of ~6.9 hours, is of comparable power to the 5.3 hour peak. It is thus possible that either 
one of these periods (i.e 5.3 or 6.9 hours) may correspond to the orbital period but it is difficult to 
discriminate between the two based solely on our data. 
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Figure 4. The binned pulse profile folded at 1559.2 seconds corresponding to the white dwarf spin period. 
The complete cycle is divided into 25 bins. The ordinate represents the relative magnitude and the abscissa 
represents the fractional phase. The spin cycle is repeated once for clarity. 



However, we propose that the 5.3 h period more likely corresponds to the orbital period of 
the binary system. First, the period of 5.3 ho ur is in reasonably good agreement with the orbital 
period of 5.54 h inferred by Schlegel (2005) from the separation of spin-orbit sideband peaks in 
the power spectrum of the X-ray data. Moreover a peak-like feature, albeit of weak power, is 
seen at a frequency 0.74478 mHz in Figure 2. If the 5.3 h period is correct, then such a feature 
is expected to arise as a result of the u + 2Q. sideband peak. On the other hand, if the 6.9h 
period is correct, then based on the same argument advanced above, we should expect a feature 
at 0.73240 mHz which is however not seen. In summary, there are indications from our data for 
a 5.3 h orbital period in the WX Pyx system. The evidence for this is not overwhelming, but 
when results from the X-ray data are also considered collectively, the 5.3 h period appears to be 
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20 




Frequency (jLtHz) 

Figure 5. The LS periodogram near the orbital frequency is shown. It may be noted that there is no 
significant power above the noise level at the first harmonic of the orbital frequency. 

a realistic possibility. We thus assume that 5.30 + 0.02 hour as the orbital period and adopt it, 
in following subsections, to determine some of the system parameters of WX Pyx. The orbital 
phase plot of WX Pyx is given in Figure 6 using a precise value of 5.297 hour(19072 second) for 
the orbital period. The time of the maximum light of the orbital modulation is found to be RJD 
54327.1897 +0.0002 d. 



3.2 Spectral type of the secondary 

It is well established that the secondary star in cataclysmic variables, with period less than 6 
hours, resemble late-type main sequence stars i.e. stars of spectral type K or M. However, the 
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Figure 6. The binned phase plot of 5.30 hour corresponding to the orbital period of WX Pyx is shown. 
The complete cycle is divided into 15 bins. The ordinate represents relative magnitude and the abscissa 
represents the fractional phase. The profile is repeated once for clarity. 



secondary stars of CVs have been shown t o be a little cooler than the isolated main sequence 
star of same mass and radius (Knigge (2006) and Figure 7 therein; hereafter K2006) An empir - 
ical relationship between orbital period and spectral type has been obtained ( B eu ermannl (1 1 99 8 ) : 
K2006). Using the empirical relation discussed in K2006, and assuming an orbital period of 5.3 
hours, we estimate the spectral type of the secondary to be constrained in the range M2 + 2. 



10 



V. H. Joshi et. al. 



3.3 Distance to the source 



Since the spectral t ype of the syst em is known, one can estimate the distance to it using the method 



first suggested by Bailey ( 1996b . This inter-relates the surface brightness at wavelength A, 



radius of the star R, distance d and observed flux f^, by following relation 



2 



Assuming that the late type secondary star is the dominating contributor to the K band flux 
and expressing flux in magnitudes we can write, 



S t = niK + 5 - 5logd + 5log I — j (2) 



where ra# is the apparent magnitude and S k is the surface brightness in the K band . The ra- 
dius R of the star can be calculated from the Roche-lobe geometry. However, the above approach 
gives only a lower limit of the distance instead of distance itself because of our assumption that the 
secondary star is the sole contributor to the K band flux. K2006 discussed the relations between 
mass-radius and spectral type-orbital period and used these in the above equation to derive the 
simpler form given below viz. 



5log(d lim ) = m K - M K {P orb ) + 5 (3) 



where dim is the lower limit to the distance, is observed K band magnitude and Mk is the 
absolute magnitude as a function of orbital period. 

From K2006 we obtain an absolute K band magnitude Mk = 5.15 for a 5.3 hour orbital 
period. The K s band magnitude of WX Pyx from 2M ASS data (K s mag . = 14.817) was converted 



to the CIT system using the transformation given by Carpenter (2001) resulting in a transformed 
CIT K magnitude of 14.857. The use of this value in the last equation gives a lower limit to 
the distance as 874 pc. However, as pointed out by K2006, the assumption that the secondary 
contributes predominantly to the K band flux is not necessarily true in a general sense. K2006 
studied the objects with well-determined distances based on the parallax method and gave an 
offset value in magnitudes which is to be applied to better estimate the distance to the object. 
Applying this offset of 1.22 magnitude for the K band we infer the distance to the source to be ~ 
1.53 kpc. 
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3.4 Mode of accretion onto the white dwarf 



In their first optical follow-up observations, Hertz et al. (1990) analyzed the optical spectrum of 



WX Pyx and calculated the equivalent width of the HyS emission line. F rom this, they estimated 



the V band absolute magnitude of the disc M„ using the relation given by [Patterson ( 19841) which 



correlates these two parameters . They compared M„ with the observed V band magnitu de to infer 
that th e distance to WX Pyx is ~ 1.6 kpc. We keep in mind that the relation given by [Patterson 



( 1984) holds for IPs with well-developed discs around the white dwarf. Interestingly the distance 
esti mate in this work of ~ 1.53 kpc, calculated using a different approach, closely matches that 
of Her tz et al. (1990). The agreement/consistency of two independent methods yielding similar 



distance estimates to WX Pyx suggests that the inherent assumption in the Hertz estimate (i.e. 
the system has a disc) is valid. Thus it would appear that the WX Pyx contains a partial disc 
around the white dwarf and it is mildly suggested that the mode of accretion is disc-fed. Further 
supporting evidence for the presence of the disc comes from the modulation at orbital period. The 
existence of such a modulation can be explained only by varying aspect of the bright spot formed 
at the location of stream-disc interaction which, again, requires the presence of the disc. 



3.5 Constraints on the angle of inclination of the system 

It is expected that CVs show a ellipsoidal variation in their near-IR light curves at the first har- 
monic of the orbital period. Such a variation would arise due to the aspect variation of the 
Roche-lobe filled distorted secondary star. Fixing the other known parameters of the system one 
can infer the angle of inclination i from the the shape of the phase curve. As seen in the Figure 1, 
the periodogram does not show significant power at the first harmonic of the orbital period shown 
as dashed line at the low frequency end of the spectrum. This suggests that the inclination angle 
of the binary system is low. 

To estimate the upper limit of the inclination angle i we produced a binned phase curve at 
the first harmonic of the orbital period. This phase curve shows a maximum scatter of +0.03 
magnitudes. In order to estimate the upp er limit of ; we produced the synthetic phase curve 



using Wilson-Devinney light curve code (Wilson & Devinnevl ll971h. Limb darken ing values 
were obtained from interpolation of data in the table provided by I Van Hamme d 1993b . The mass 



ratio q was varied between 0.2 and 0.6 and the appropriate inclination angle was obtained for 
each value of q such that the synthetic light curve shows a similar scatter as that observed in the 
data. We find that for any value of q in the chosen range, the inclination angle i is less then 30°. 
This is suggestive that the inclination angle of the system is small. 



3.6 Mass transfer rate 



The mass transfer rate M of the system can be estimated by comparing the accretion luminosity 
L acc and the released gravitational potential energy. When the mass is transferred from the inner 
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Lagrangian point to the white dwarf surface, the gravitational potential energy is released and 
radiates away, largely at X-ray wavelengths in the case of magnetic CVs. 



To estimate the X-ray luminosity of WX Pyx we used the flux values given by Evans & H ellier 



(2007). These authors modeled the X-ray spectrum of WX Pyx using a combination of a non- 
thermal hard component and a thermal soft component. The unabsorbed bolometric X-ray flux 
corresponding to the soft and hard components are estimated to be 6.0 x 10 -13 ergs s _1 cirT 2 and 
7.51 x 10~ 12 ergs s cirT 2 respectively giving a total X-ray flux of 8.11 x 10~ 12 ergs s _I cirT 2 . 
This flux is converted into the X-ray luminosity by the relation L x = 2nd 2 F x where L x is the X- 
ray luminosity in ergs s , F x is the X-ray bolometric flux in ergs s _I cirT 2 and d is the distance. 
Here, the factor 2n is used instead of An to take into account that the X-ray emitting shock region 
is close to the white dwarf surface and half of the total emission is blocked by the stellar surface. 
Using a value of d — 1.530 kpc gives the X-ray luminosity L x = 1.1 X 10 33 ergs s _1 . Estimation 
of L acc is a complex process requiring information about various aspects including the IR to UV 
emission from the accreting material, absor ption effect a t UV and s oft X-ray wavel engths and the 
contribution from X-rays above 10 keV (lWarnerlll995l) . However I Warner] (11995b indicates that 
L acc ma y be estimated, within an uncertainty of a factor of 2, by multiplying L x by 50. Therefore, 
L x was multiplied by this factor to estimate L acc which in turn is converted into the mass transfer 
rate M using L acc = GMwdM/Rwd where L acc is the total accretion luminosity and Mwd and 
Rwd are the mass and radius o f the white dwarf respe ctively. Mwd was chosen to lie in the range 
1.31 to 1.4 M as inferred by Evans & Hellierl (12007 ) from modeling of the hard component of 
the X-ra y data. Rwd was calculated from the empirical mass radius relationship for white dwarfs 
given by Prin gle & Webbinkl (119751) . We obtain the mass transfer rate to lie in the range 0.6 to 
l.OxlO 17 g s -1 or 0.95 to 1.6xl0~ 9 M o yr _1 . This v alue of M is fo und to be fairly typical for IPs 
as has been compiled for several other such objects dWarnerlll995l) . 



3.7 Magnetic moment of the white dwarf 



Assuming a magnetic white dwarf with spherically symmetric accretion, the magnetospheric 
radius r mag is defined to be the radius at which the magnetic p ressure balances the ram pres sure 
of the infalling material. Following the calculation given in [Frank. King & Ra ine (2002) the 
magnetospheric radius is given by 



2ttV 



2, ,4 



H 2 GM WD M 2 



1/7 



(4) 



where yU and //o are magnetic moment of white dwarf and permeability of free space respect- 
ively. However, if the white dwarf is accreting via a disc, the true magnetospheric radius R mag 



can be estimated from R mag — 0.5 r mag . 



If the white dwarf of the IP is in spin equilibrium it is assumed that the magnetospheric radius 
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Rmag is very close to the co-rotation radius R co (e.g. Norton et al 2004). Here R co is the radius at 
which the accreting material in local Keplerian motion co-rotates with the magnetic field of the 
white dwarf. Therefore, 



R 



mag 



Rc 



(GM WD P 
An 2 



1/3 



(5) 



Comparing equations (4) and (5), we get 



(6) 



For WX Pyx, /i is thus estimated in the range 1 .9 to 2.4x 10 33 G cm 3 for a Mwd value ranging 
between 1.31 an d 1.4Mp. This result is consistent with the magnetic moments of several other 
IPs estimated by I Warner! d 1995b . 



Nort on et al.1 Q2008) have simulated the accretion flow of magnetic CVs in a detailed study. 



They have shown that for a fixed orbital period and mass ratio, different regions in a diagnostic 
plot of the spin period versus magnetic moment, correspo nd to differen t geo metries for the ac- 
cretion viz. disc, stream, ring and propeller (Fig 1 and 2 in [Norton et al. ( 2008 )). For the orbital 
period and mass ratio specific for WX Pyx, we find that the estimated value of magnetic moment 
along with the spin period falls well inside the region corresponding to a disc type accretion. This 
also supports disc-fed accretion as the favored mode of accretion in WX Pyx. 



4. Summary 



We have presented near infrared 1 .25 /vm / band photometry of WX Pyx covering a total duration 
of 22 hours spanning six nights from December, 2007 to February, 2010. Our motivation was to 
determine the spin period of the object from IR observations, an exercise which has not been done 
earlier, and if possible to also determine the orbital period. The frequency analysis of the light 
curve clearly indicates the presence of a spin period of 1559.2 ± 0.2 seconds for the white dwarf. 
However, the orbital period is less robustly determined. From the IR observations, a likely peak 
at 5.30 + 0.02 hour is seen in the power spectrum of the object which is argued to be the orbital 
period of the system. Subsequently, estimates are made of some of the physical properties and 
parameters of the system viz. the spectral type of the secondary star, the distance to the object, the 
mode of accretion, the angle of inclination of the system, the mass transfer rate and the magnetic 
moment of the white dwarf. 
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